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ABSTRACT 



Q 

Context. Accurate physical properties of eclipsing stars provide important constraints on models of stellar structure and 
evolution, especially when combined with spectroscopic information on their chemical composition. Empirical calibrations of 
the data also lead to accurate mass and radius estimates for exoplanet host stars. Finally, accurate data for unusual stellar 
^ , subtypes, such as Am stars, also help to unravel the cause(s) of their peculiarities. 

$_( ■ Aims. We aim to determine the masses, radii, effective temperatures, detailed chemical composition and rotational speeds for 
jy^ the Am-type eclipsing binaries SW CMa ( A4-5m) and HW CMa ( A6m) and compare them with similar normal stars. 

Methods. Accurate radial velocities from the Digital Speedometers of the Harvard- Smithsonian Center for Astrophysics were 
combined with previously published uvby photometry to determine precise physical parameters for the four stars. A detailed 
I ^ abundance analysis was performed from high-resolution spectra obtained with the Nordic Optical Telescope (La Palma). 
^ . Results. We find the masses of the (relatively evolved) stars in SWCMa to be 2.10 and 2.24 Mq, with radii of 2.50 and 3.01 Rq, 
' while the (essentially zero-age) stars in HWCMa have masses of 1.72 and 1.78 Mq, radii of 1.64 and 1.66 -Rq - all with errors 
, well below 2%. Detailed atmospheric abundances for one or both components were determined for 14 elements in SWCMa 
G^ ■ ([Fe/H] = -I-0.49/-I-0.61 dex) and 16 in HWCMa ([Fe/H] = -I-0.33/-I-0.32 dex); both abundance patterns are characteristic of 
' metallic- line stars. Both systems are well fit by current stellar evolution models for assumed bulk abundances of [Fe/H] = 
' -1-0.05 and -1-0.23, respectively ([a/Fe] = 0.0), and ages of ~700 Myr and 160 Myr. 

T-H ■ Key words. Stars: binaries: eclipsing - Stars: fundamental parameters - Stars: abundances - Stars: chemically peculiar - Stars: 
■ individual: SW CMa - Stars: individual: HW CMa 

^ 1. Introduction of main-sequence A stars are known to be metallic-lined 

■ - - - , ... , . . ■ , , 1 .1 (Am). The Am stars are found overwhelmingly to be in 
Among- the eclipsmg bmaries with the best determmed , . , , ai .nrrTmi rm i ^ • ju 
_ , , . . , . bmary systems (e.g., Abt|1961J. 1 hey are characterized by 
lundamental properties approximately two dozen are oi ,. , , . , ,. , r. • n- i , 
. ; , , , „„rw peculiar abundance patterns including a deficiency ot hght 
spectral type A (see Torres et al.i2010i. Roughly 20-30% , , , c jr. r. ^ f i i 
^ ' ° elements such as Sc and Ca, enhancement oi iron-peak ei- 

Send offprint requests to: G. Torres, ements, and a strong overabundance of Sr and some rare 

e-mail: gtorres@cfa.harvard.edu earths, sometimes by an order of magnitude or more. The 

* Based on observations carried out with the Nordic Optical enhancement generally increases with increasing atomic 

Telescope (NOT) at La Palma, the 50 cm Stromgren Automatic number. Their typically slower axial rotation and stronger 

Telescope (SAT) at ESO, La SiUa, the 1.5 m Wyeth reflec- nietal lines compared to normal A stars tend to facilitate 
tor at the Oak Ridge Observatory, Harvard, Massachusetts, radial-velocity measurements necessary for accurate 

USA, and the 1.5 m TiUinghast reflector at the F. L. Whipple ^^^^ determinations. It is therefore not surprising that 
Observatory, Mt. Hopkins, Arizona, USA. , -,1 a , , ^ ,i 

SZ—, .r-r-T^ , ^ , systems with Am stars are over-represented among the 

tables I A. II and I A. 21 as well as Appendices B and C will be . r n i ir r i a 

., , , . , ^ .J. ^ ^, r^rsc^ ■ ri i bcst-measurcd eclipsing binaries: lully halt of the A-type 

available m electronic form at the CDS via anonymous ftp to jr- o j j 

130.79.128.5 or via .http://cdsweb.u-strasbg.fr/Abstract.html, systems for which we know the absolute masses and radii 
*** Deceased. 
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with relative errors less than 3% contain at least one com- 
ponent showing these abundance abnormalities. 

Although significant progress in understanding and 
modeling these chemically peculiar objects has been made 
in recent years, important questions regarding the role of 
rotation and the origin of the detailed patterns of element 
enhancement and depletion are not yet completely set- 
tled (see, e.g., Abt 1^0001 Bohm-Vitense [2003 Vick et al. 
I2010p . Further observational constraints are highly desir- 
able, particularly ones in which very accurate mass and ra- 
dius determinations are accompanied by a detailed abun- 
dance study of the components. 

Beyond a qualitative assessment of their metallic-lined 
nature, individual element abundances are known for pre- 
cious few of the dozen or so Am binaries with well deter- 
mined properties (see, e.g., Lyubimkov et al. I1996|) . Thus, 
little can be said about whether or how the detailed abun- 
dance patterns that seem to vary to some degree from 
star to star depend on global properties such as mass, ef- 
fective temperature, or surface gravity. In this paper we 
present extensive new observations of two eclipsing Am 
systems, SW CMa and HW CMa, with the goal of deter- 
mining not only accurate absolute dimensions but also 
individual abundances for all four stars. It is hoped that 
these results and others like them will present an oppor- 
tunity for advancement on some of the issues mentioned 
above. Figure [1] places the two systems in the context of all 
other main sequence eclipsing binaries in the same mass 
regime that have the best-determined masses and radii. 
Am stars are marked with open circles. SWCMa is seen 
to be among the more evolved cases, while HWCMa is 
the least evolved in this mass range. 

We describe below our spectroscopic observations of 
these systems, and an analysis of those data along with 
previously-published high-precision four-color Stromgren 
photometry in order to obtain accurate masses and radii. 
For SWCMa our determinations substantially improve 
upon previous estimates; for HW CMa they represent the 
first such measurements. We conclude with a detailed com- 
parison with models of stellar and tidal evolution. 

2. SWCMa 

The observational history of SWCMa (HD 54520, 
HIP 34431, V = 9.15, A4-A5) has been summarized by 
Clausen et al. (|2008p . Since its discovery as an eclipsing bi- 
nary by Hoflmeister (jl932p , relatively little progress in un- 
derstanding its physical properties was made in more than 
six decades until the most recent study by Lacy (jl997p . In 
that work photoelectric UBV light curves were combined 
with new spectroscopy to yield the first precise determi- 
nation of the masses and radii for the components. The 
orbit is eccentric, and apsidal motion was detected by Lacy 
(|1997p and subsequently refined by Clausen et al. 12008), 
although the apsidal period is very long and still poorly 
constrained {U = 14 900 ± 4700 yr). 



For the analysis in this paper we adopt the linear 
ephemeris derived by Clausen et al. (|2008p . 

Min I (HJD) = 2 446 829.6482(1)-!- 10^091988(5) x E , 

in which the uncertainties in the period and reference 
epoch are indicated in parentheses in units of the last deci- 
mal place. This corresponds to the eclipse of the more mas- 
sive and slightly hotter component. The secondary eclipse 
occurs approximately at phase 0.31. 

2.1. Radial velocities and spectroscopic orbit 

The spectroscopic observations of SW CMa used here for 
the radial-velocity determinations were carried out at the 
Harvard-Smithsonian Center for Astrophysics (CfA) be- 
tween 1988 April and 1989 March. A total of 39 spec- 
tra were gathered with the 1.5 m Wyeth refiector at 
the Oak Ridge Observatory (Harvard, Massachusetts). 
Two additional spectra were taken using the 1.5 m 
Tillinghast reflector at the F. L. Whipple Observatory 
on Mt. Hopkins, Arizona. Nearly identical spectrographs 
("Digital Speedometers" (DS); Latham [1992p were used 
on each telescope, equipped with photon-counting Reticon 
detectors that recorded a single echelle order 45 A wide 
centered at a wavelength near 5187 A (Mgib triplet). The 
resolving power of these instruments is A/AA w 35 000, 
and the signal-to- noise (S/N) ratios of the 41 spectra we 
collected range from 10 to 25 per resolution element of 
8.5 kms-i. 

Radial velocities were obtained with TODCOR, a 
two-dimensional cross-correlation technique introduced 
by Zucker & Mazeh (|1994p . Templates for the cross- 
correlations were selected from a library of synthetic spec- 
tra based on model atmospheres by R. L. Kurucz, which 
has been described by Nordstrom et al. (|1994p and Latham 
et al. (|2002p . The parameters of these templates are the ef- 
fective temperature {Tes), metallicity ([m/H]), rotational 
velocity (wsinz, when seen in projection), and surface 
gravity (log 5). 

The optimal templates for the two components 
were determined by running extensive grids of cross- 
correlations and selecting the combination of parameters 
yielding the highest correlation averaged over all expo- 
sures, weighted by the strength of each spectrum (see 
Torres et al. I2002p . The surface gravities were held fixed 
at \ogg = 4.0 for both stars, close to the final values in 
Sect. 12.41 Temperatures and rotational velocities were op- 
timized for a range of different (fixed) metallicities, and 
the first hint that the stars have anomalously strong metal 
lines came from the fact that the highest cross-correlation 
value was obtained for the largest metallicity available in 
our grid, [m/H] = +0.5, for both components. It is possi- 
ble that a better match to the observed spectra could be 
found for an even higher metallicity, and from the trend of 
the improvement with [m/H], the effect appears somewhat 
more pronounced for the secondary. 

By interpolation the optimal temperatures were found 
to be 8210 K and 7860 K, and the best t^sini values are 
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29 kms^^ and 15 kms~^. However, other photometric and 
spectroscopic estimates discussed later point toward very 
similar Tcs values for the two stars, as opposed to the 
350 K difference implied above. Our primary star Tes es- 
timate is in agreement with those determinations, but the 
secondary value seems too cool. Due to the limited wave- 
length coverage of the CfA spectra, metallicity is strongly 
correlated with effective temperature so that an almost 
equally good fit to the spectra can be found for different 
template combinations in which the stronger lines corre- 
sponding to a metal-rich composition can be made to ap- 
pear weaker by increasing the temperature, and vice versa. 
Therefore, even though we cannot go any higher in [m/H] 
in our optimization scheme, a satisfactory match to the 
observed spectra can still be found by adjusting the tem- 
perature. This is likely the explanation for the excessively 
cool Tes we derive for the secondary from these spectra, 
which supports the notion that this star may have a more 
anomalous (generally more metal-rich) composition than 
the primary. Experience shows that this degeneracy has 
little impact on the radial velocities themselves, although 
it does of course limit our ability to simultaneously infer 
accurate values of [m/H] and Toff. For the velocity deter- 
minations we adopted [m/H] = -1-0.5 for both stars, along 
with temperatures of 8250 K and 7750 K and vsini pa- 
rameters of 30 kms~^ and 16 kms~^ for the primary and 
secondary. These are the values in our grid nearest to the 
optimal estimates mentioned earlier. We do not attach any 
astrophysical meaning to these T^ff values, but treat them 
merely as adjustable parameters intended to provide the 
best radial velocities, as free as possible from biases. 

The stability of the velocity zero point during the pe- 
riod of observation was monitored by taking exposures 
of the dusk and dawn sky. Small run-to-run corrections 
were applied in the manner described by Latham (|199 2!). 
Additional systematic errors may occur because of the 
narrow wavelength coverage of the CfA spectra and the 
fact that lines move in and out of the spectral window 
as a function of orbital phase and the heliocentric correc- 
tion (see Latham et al. I1996|) . To investigate this we per- 
formed experiments with simulated double-lined spectra 
following the procedure detailed by Torres et al. (|1997p . 
and corrections based on them were then applied to the 
raw velocities. For SW CMa these corrections are typically 
small (< 1 kms~^), and their net effect is to increase the 
absolute masses by slightly less than 0.4%, and the radii 
by about 0.1%. They are shown in Fig. [21 where the sys- 
tematic pattern is obvious, particularly in the lower panel. 

The final radial velocities referred to the heliocentric 
frame are presented in Table lA.ll of Appendix [2] (avail- 
able electronically), and include all corrections mentioned 
above. The elements of the orbital solution we obtain from 
them are listed in Table [1] along with derived quantities 
including the minimum masses and semimajor axes. The 
measurements and the fitted orbit are shown in Fig. [3] to- 
gether with the residuals. The larger rms residual for the 
primary, despite being the brighter star, is explained by 
the significantly higher rotational broadening compared to 



Table 1. Spectroscopic orbital solution for SWCMa. 



Parameter 



Value 



Adjusted quantities: 
Kp (kms-i) 
Ks (kms^^) 
7 (kms^^) 

Adopted quantities: 
P (days) 

Ti (HJD-2 400 000) ^ 

Derived quantities: 
Mp sin^ i (Mq) 
MsSin^i (Mq) 
q = M^/Mp 
ap sin i (10^ km) 
as sini (10^ km) 
asini [Rq) 

Other quantities pertainini 

Nobs 

Time span (days) 
Up (kms~^) 
as (kms""*^) 



82.09 ±0.35 
87.37 ±0.22 
+42.57 ±0.17 
0.3174 ±0.0011 
163.26 ±0.51 

10.091988 
46829.6482 

2.237 ±0.014 
2.102 ±0.018 
0.9396 ± 0.0048 
10.803 ± 0.046 
11.498 ±0.029 
32.057 ± 0.076 

I to the fit: 

41 
333 
1.93 
1.19 



Longitude of periastron for the more massive star. 
^ Time of central primary eclipse (eclipse of the more mas- 



sive star). 



the secondary. We note that while our secondary velocity 
semi-amplitude agrees with the value reported by Lacy 
(|1997p . our primary semi-amplitude is slightly larger (by 
2%, a 1.6-cr difference). 

In addition to the radial velocities, we used TODCOR 
to derive the light ratio at the mean wavelength of our 
observations (5187 A), following Zucker & Mazeh (|1994|) . 
We obtain Ls/ Lp = 0.68 ± 0.04 (secondary/primary). 

2.2. Chemical abundances 

In order to obtain reliable estimates of the spectro- 
scopic properties, a spectrum of SWCMa was obtained 
on UT 2008 November 17 (mid exposure at HJD 
2 454787.65730, orbital phase 0.5472) using the FIES in- 
strument (Frandsen & Lindberg ll999p on the 2.5 m Nordic 
Optical Telescope (NOT) on La Palma, with a resolv- 
ing power of A/ A A w 46 000. The signal-to-noise ratio 
achieved is about 60 per pixel in the 6160 A region, and 
the wavelength coverage is approximately 3600-7400 A 
recorded in 78 echelle orders. Reductions used the IRAF- 
based FIEStool packag43 with nightly flatfield, bias, and 
dark frames, as well as thorium-argon exposures taken im- 
mediately before and after the science exposure. 



^ See |http : / /www . not . iac . es] for details on FIES and 
FIEStool. 



4 



G. Torres et al.: SWCMa and HWCMa 



The spectroscopic analysis was carried out with the 
IDL-based VWA analysis tool, extended to the analysis of 
double-lined spectra (for details of the procedure, we refer 
the reader to Bruntt et al. 2004, 2008, 2009). Briefly, VWA 
uses the SYNTH software (Valenti & Piskunov HM^ to 
generate synthetic spectra. Atmosphere models are inter- 
polated from the recent grid of MARCS model atmo- 
spheres (Gustafsson et al. I2008p . which adopt the solar 
composition by Grevesse et al. (|2007p . Line information 
was taken from the Vienna Atomic Line Database (VALD; 
Kupka et al. I1999p . In deriving abundances relative to the 
Sun, the log gf values were adjusted in such a way that 
each measured line in the Wallace et al. (|1998p solar atlas 
reproduces the atmospheric abundances by Grevesse et al. 
(PPUTI) . 

The number of unblended lines for the analysis is some- 
what limited in this case due to the relatively large ro- 
tational broadening. We restricted the measurements to 
lines with equivalent widths smaller than about 65 mA in 
the composite spectrum (corresponding to about 110 mA 
for the primary and 160mA for the secondary), and the 
surface gravities were held fixed to preliminary values near 
those given in later Table HI We assumed a light ratio of 
Ls/Lp = 0.69 {V band), from the light-curve solutions 
described below. The effective temperatures and micro- 
turbulent velocities (umic = 2.70 kms~^ and 2.85 kms~^ 
for the primary and secondary, with estimated errors of 
0.8 kms~^) were adjusted to yield consistent values for 
the abundances of Fei and Fell, and to avoid correlations 
with the equivalent width and excitation potential. The 
results for T^s indicate similar values for the primary and 
secondary of 8200 ± 150 K and 8100 ± 150 K, respectively, 
corresponding approximately to spectral types A4 or A5. 

Individual abundances on the scale of the Grevesse et 
al. (|2007[) solar abundances were determined separately 
in one or both stars for 14 species. These results are 
listed in Table [H along with the number of spectral lines 
used in each case. Abundances based on fewer than three 
lines are somewhat more uncertain. NLTE effects become 
significant in stars with temperatures hotter than about 
8000 K, so appropriate corrections to the abundances of 
Fel (-1-0.13 dex and -1-0.11 dex for the primary and sec- 
ondary) have been included in the values we report, fol- 
lowing Rentzsch-Holm (1996J). The average iron abun- 
dances for the components are [Fe/H] = -1-0.49 ± 0.15 and 
4-0.61 ± 0.15, where the uncertainties include a contri- 
bution from the errors in temperature and Umic- Figure 2] 
displays the abundance pattern for both stars in SW CMa. 
Although some of these determinations are rather uncer- 
tain due to the small number of lines and limited signal- 
to-noise ratio of our spectroscopic material, the enhanced 
iron-peak abundances compared to calcium, and especially 
the strong overabundance of heavy elements such as Sr, 
Y, and Ba, are typical of Am stars and suggest that both 
components are chemically peculiar. The iron abundance, 
as well as the Fe/Ca difference (0.44 and 0.81 dex), is 
somewhat larger for the secondary than the primary, con- 
sistent with the hints we saw in the CfA spectra. The pro- 



Table 2. Abundances for the components of SWCMa. 



Primary 



Secondary 



Element 


[X/H] 


a 


iV 


[X/H] 


a 


iV 


Nal 


+0.57 




1 








Mgl 


+0.21 




1 








Si I 








-0.51 




2 


Sill 


-f0.20 




2 








Cai 


-f0.05 


0.17 


3 


-0.20 


0.16 


4 


Tin 


+0.02 


0.21 


3 


+0.37 


0.12 


3 


Cm 


+0.15 




2 


+0.19 




1 


Fei 


+0.50 


0.09 


17 


+0.59 


0.10 


22 


Fell 


+0.47 


0.09 


11 


+0.63 


0.11 


15 


Nil 


+0.51 




2 


+0.42 




2 


Zni 


+1.08 




2 


+1.06 




1 


Srii 


+1.00 




1 


+1.29 




1 


Yii 


+1.02 




1 


+1.11 




1 


Ban 


+1.40 


0.21 


3 


+1.54 




2 



Note: These abundances are based on Vmic = 2.70 kms~^ 
and 2.85 kms"'^ for the primary and secondary, respectively, 
and [Fe/H] = +0.45 for the model atmospheres applied in the 
analysis. Fe I abundances include corrections for NLTE. A'^ rep- 
resents the number of lines of each element, and a is the scatter 
of the measurements (listed only if A'^ > 2). 



jected rotational velocities were measured from the FIES 
spectrum by synthesizing line profiles for a dozen isolated 
lines and comparing them with the observed spectrum, 
seeking to minimize the residual differences. The adopted 
macroturbulence velocity was ^rt — 8kms^^ for both 
components, extrapolated from the trend with tempera- 
ture reported by Gray (|2005p . We obtained wsini values 
of 24.0 ± 1.5 km s^^ and 10.0 ± 1.0 km s^^ for the primary 
and secondary, respectively. That these rotational veloci- 
ties are much slower than typical for A stars in the field 
is a general characteristic found in other Am stars as well 
(see, e^., Abt & Morrell ITM51 Abt [MM Fossati et al. 



2.3. Photometric elements 

For this analysis we make use of the differential photom- 
etry for SW CMa in the Stromgren system that has been 
reported previously by Clausen et al. (|2008|1 . As described 
there, the uvby light curves consist of 820 observations in 
each band collected with the 0.5 m Stromgren Automatic 
Telescope at ESQ (La Silla), on 89 nights during five peri- 
ods between 1987 February and 1991 March. The average 
uncertainty per differential observation is about 5 mmag 
in vby and 6 mmag in u, and the measurements scatter 
at these same levels throughout all phases, indicating that 
the components of SW CMa are constant within the pre- 



Note, however, that slow rotation in an A star is not al- 
ways associated with chemical pecuharity (see, e.g., Fekel et 
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Table 3. Adopted photometric elements for SWCMa. 



i (°) 88.59 ±0.20 

e cos cj -0.30313 ± 0.00009 

esinuj +0.0883 ± 0.0017 

e 0.3157 ± 0.0004 

Lu C) 163.76 ±0.29 

rp + T-s 0.1718 ±0.0010 

k 0.827 ±0.018 

r-p 0.0940 ± 0.0006 

rs 0.0778 ± 0.0013 



y h V u 

Js 1.013 1.012 0.996 1.020 

±3 ±2 ±3 ±4 

Ls/Lp 0.693 0.692 0.681 0.697 

±30 ±18 ±15 ±33 



Note: The individual flux and luminosity ratios are based on 
the mean stellar and orbital parameters. 

cision of the data. Compared to the UBV light curves in 
the earlier study by Lacy (|1997p . the present observations 
are roughly twice as numerous and provide better orbital 
coverage, particularly at the shoulders of the minima. 

Light-curve solutions were performed with the 
JKTEBOP code (Southworth et al. I2004al [2004b| . which 
is an updated and expanded version of the original EBOP 
program based on the Nelson-Davis-Etzel model (Nelson 
& Davis mill Etzel UMT]) . This model represents the bi- 
nary components as biaxial ellipsoids and applies a sim- 
ple bolometric reflection prescription, but is perfectly ad- 
equate for the stars in SW CMa given their very small de- 
formation (oblateness of 0.0012 for the primary and 0.0008 
for the secondary, both much smaller than the tolerance 
for this model of 0.04; see Popper & Etzel .198,1) . Each of 
the uvby light curves was analyzed independently adopt- 
ing the ephemeris given in Sect. [5J with equal weight as- 
signed to all observations. The main parameters solved 
for are the radius ratio between the secondary and the 
primary (k = Vs/rp), the sum of the radii (rp + r^), 
the inclination angle (i), the geometric factors esinw and 
ecosw, the central surface brightness ratio of the sec- 
ondary in terms of the primary (Js), a photometric scale 
factor (the magnitude at quadrature), and a correction 
to the phase of the primary minimum. The mass ratio 
was held fixed at the value determined spectroscopically 
(Sect. [2TT|) . Gravity darkening coefficients, yp and were 
computed from a simple Planck approximation (see, e.g., 
Martvnov I1973p . Initial solutions were carried out with 
limb-darkening (LD) coefficients interpolated from theo- 
retical calculations (van Hamme 119931 and ATLAS ver- 
sion of the Claret 120001 tables) according to the effective 
temperatures and surface gravities for stars as determined 
hereH 

^ Solar metallicity was adopted in interpolating all limb- 
darkening coefficients for consistency in the comparison be- 



We found good agreement between the results from 
the uvby passbands, although small systematic differences 
were noticed between fits using LD coefficients from differ- 
ent sources, as described in detail in AppendixlB] (available 
electronically). In view of this, and out of concern that im- 
posing LD coefficients from theory might introduce subtle 
biases in the geometric elements, for the final fits we chose 
to leave the LD coefiicients free (linear law) , with the only 
constraint that they be equal for the primary and sec- 
ondary given that the temperatures are also very nearly 
the same. Additional experiments with LD free were car- 
ried out to explore the effect of including third light (^3) as 
an extra parameter. For all four passbands we found that 
£3 was negligible compared to its error, and we conclude 
third light is not significant. 

The results are presented in Table |3l in which the 
values from the separate passbands have been averaged 
together, with weights inversely proportional to the rms 
residual of each solution. The uncertainties listed in the 
table were obtained by running 1000 Monte Carlo simula- 
tions with JKTEBOP, and we consider them to be more 
realistic than the formal (internal) errors from the solu- 
tions presented in Appendix [BJ For the relative radii we 
also computed uncertainties by propagating the errors in 
Tp + Ts and k following Torres et al. (|2000p . This procedure 
assumes the sum and ratio of the radii are uncorrelated, 
and results in uncertainties that are more nearly equal 
(0.0011 and 0.0010 for Vp and r^, respectively) than those 
from Monte Carlo. However, the Monte Carlo simulations 
indicate that in this case Vp + Vs and k are in fact cor- 
related, invalidating the error propagation approach for 
SW CMa. The radiative quantities in Table [3] were com- 
puted from additional solutions in which the geometry was 
held fixed to the weighted average from the four bands. 
The eccentricity and longitude of periastron are more pre- 
cise than, but consistent with those obtained spectroscop- 
ically, differing by 1.5a and 0.9cr, respectively. The light 
ratios in b and y (Ls/Lp = 0.69 ± 0.02 and 0.69 ± 0.03) 
are in excellent agreement with the value from our spec- 
troscopic observations in Sect. [Hi] (0.68 ± 0.04 at 5187 A). 
According to the final light-curve elements, 72% of the 
light of the primary is blocked at the deeper minimum. 
The secondary eclipse is nearly total: 97% of the y-band 
light of that star lost at that phase. 

The light curves along with our best fit models appear 
in Fig. [51 and expanded views of the primary and sec- 
ondary minima are shown in Fig. [S] and Fig. [7] The O — C 
residuals are displayed in the bottom panels. 



tween the different LD tables, given that the van Hamme 
(il993) calculations do not reach metallicities as high as those 
measured for SWCMa. However, inspection of the Claret 
1)2000 1) calculations shows that the coefficients for [Fe/H] = 
0.0 and [Fe/H] = ±0.5 differ very little. 
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Table 4. Astrophysical properties of the SW CMa system. 





Primary 


Secondary 


Absolute dimensions: 






M (Mq) 


2.239 ±0.014 


z.iU4 It U.Ulo 




3.014 ±0.020 


AQK -I- n 0/10 


logs (cgs) 


3.8298 ± 0.0065 


o.yo ( ± u.uio 


V sini (kms^^) 


24.0 ± 1.5 


10.0 ± 1.0 


Vsync (kms"^) ^ 


15.1 ±0.1 


12.5 ±0.2 


Wpsync (kmS~^) 


24.5 ±0.2 


20.3 ±0.3 


Upcri (kms"^) 


30.6 ±0.2 


25.3 ± 0.4 


Radiative and other properties: 




Tcff (K) 


8200 ± 150 


8100 ± 150 


A^bol (mag) 


0.817 ±0.081 


1.281 ±0.088 


log L/Lq 


1.566 ±0.032 


1.380 ± 0.035 


BCv (mag) 


±0.02 ±0.10 


±0.02 ±0.10 


My (mag) 


0.79 ±0.12 


1.26 ±0.13 


Distance (pc) 


576 ± 27 


V - Mv (mag) 


8.80 ±0.10 


[Fe/H] ^ 


±0.49 ±0.15 


±0.61 ±0.15 


Photometric indices: 








9.721 ±0.021 


10.119 ±0.029 




0.093 ± 0.012 


0.095 ±0.014 


mi * 


0.203 ±0.012 


0.220 ±0.013 




1.031 ±0.015 


0.987 ± 0.020 


E{b - y) 


0.023 ±0.010 



^ Measured projected rotational velocity. 
^ vsini expected for synchronous rotation. 

vsini expected for pseudo-synchronous rotation. 

Expected v sin i if synchronized at periastron. 

Values representative of the surface layers only. 
' Not corrected for interstellar absorption/reddening. 

Note: The bolometric corrections (BCv) are adopted from 
Flower (|19IIB|l . along with with T^^q = 5780 K, and A'/^^jq = 
4.73 (see Torres I2010|l . An additional error contribution of 
0.10 mag is added in quadrature to the BCv uncertainty prop- 
agated from the temperature uncertainties. 



2.4. Absolute dimensions 

Our spectroscopy and differential photometry for 
SW CMa yield masses with uncertainties of 0.6% and 0.9% 
for the primary and secondary, and radii that are good to 
0.7% and 1.7%, respectively. The values of M and 7? are 
consistent with the analysis by Lacy (|1997p . but our er- 
rors are generally smaller (by about a factor of three in 
the case of the masses) We list the masses and radii in 
Table m along with other derived properties. 

Interstellar reddening may be estimated from the 
uvbyP indices out of eclipse reported by Clausen et al. 
(|2008p . which agree well with those of Wolf & Kern 
([rggg]) . Using the calibration by Crawford (|W51) we ob- 
tain E{h — y) — 0.023, to which we assign a conservative 
uncertainty of 0.01 mag. We adopt this value in the fol- 



* We note that the uncertainty of 0.5 Mq for the primary 
mass reported in Lacy's Table 7 is likely a misprint; we have 
assumed it should read 0.05 Mq. 



lowing. The corresponding visual extinction is A{y) — 
0.10 ±0.04. Other estimates of the extinction in the direc- 
tion of SW CMa may be inferred from dust maps (which 
measure the total extinction), after proper correction for 
the distance to the system by iterations. However, we find 
those estimates to be larger and rather discrepant: from 
Hakkila et al. (jl997p we obtain A{V) — 0.16, correspond- 
ing to E{h — y) — 0.037, while the maps of Schlegel et al. 
P^gg)) suggest A{V) = 0.32, or E{h - y) = 0.077, at the 
distance to the system. 

The combined uvby indices from Clausen et al. (|2008p 
and the light ratios from our photometric solutions yield 
the individual indices for the components, also listed in 
Table 21 The uncertainties include the contribution from 
all observational errors. The individual bolometric lumi- 
nosities follow from the spectroscopic temperatures and 
the radii, and the distance was obtained from the lumi- 
nosities, the de-reddened V magnitude, and bolometric 
corrections from Flower (jl996|) . We infer a distance to the 
system of 576±27pc, and separate calculations for the pri- 
mary and secondary agree closely (giving 583 ± 36 pc and 
566 ± 37 pc, respectively), which is an indication of good 
internal consistency. Knowledge of the individual photo- 
metric indices for the components provides an opportunity 
to derive temperature estimates as a check on the spec- 
troscopic values, to the extent that existing temperature 
calibrations for normal A-type stars can be applied to ob- 
jects displaying the Am phenomenon. This appears to be 
the case, according to Netopil et al. p008p . After correct- 
ing the individual indices for the effects of interstellar red- 
dening, application of the calibration of Napiwotzki et al. 
([iggg)) leads to estimates of 8188±157K and 8168±173K 
for the primary and secondary, in good accord with the 
spectroscopic values. A 0.02 mag error in 6 — y changes 
these results by about 200 K. Using the central surface 
brightness ratio Jg in the y band obtained from our light 
curve solutions with the flux calibration by Popper (|1980p 
yields an estimate of the temperature difference between 
the stars of ATes = — 32K, suggesting the secondary 
(less massive star) is marginally hotter than the primary, 
though well within the uncertainties of our other determi- 
nations. 

The wsini values discussed in Sect. 12.11 that provide 
the best fit to the CfA spectra are 29 ± 3 kms~^ and 
15 ± 2 kms""'^ for the primary and secondary, respec- 
tively. Strictly speaking, however, these measures rep- 
resent the total line broadening, and not just rotation, 
and indeed they are systematically larger than those ob- 
tained from the FIES spectrum (24.0 ± 1.5 kms~^ and 
10.0 ± 1.0 kms""'^). This is most likely due to the fact that 
the synthetic templates used to analyze the CfA spectra 
have been calculated with a radial-tangential macrotur- 
bulent velocity appropriate for solar-type stars (Crt = 
1.5 kms~^), whereas the values of Crt expected for A- 
type stars (at least normal ones) are significantly higher 
(see, e.g., Grav 12005)) . Indeed, for the FIES analysis we 
adopted Crt = 8 kms~^. We therefore rely on the FIES 
determination of vsini, which we report in Table SI Also 
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listed there are the v sin i values predicted if the stars had 
their spins synchronized with the mean orbital motion 
(i^sync)j with the motion at periastron (wpcri), and with an 
intermediate equilibrium rate of 'pseudo-synchronization' 
(see Hut I198ip . The comparison would suggest the pri- 
mary component of SWCMa is rotating at the pseudo- 
synchronous rate and the secondary is spinning consider- 
ably more slowly than pseudo-synchronous. We note here 
also that Lacy (|1997p measured rather different values for 

V sin I of 30±2 kms~^ for the primary and 21±3 kms~^ for 
the secondary, which would lead to different conclusions 
regarding synchronization. 

The Am nature of both stars is supported by the abun- 
dance pattern we see, most notably that the heavy ele- 
ments are overabundant by an order of magnitude or more 
compared to solar. Enhancements seem rather similar for 
the two components, with a few exceptions noted earlier. 

3. HWCMa 

The long-period eclipsing binary HWCMa (HD 54549, 

V = 9.19, A6, P = 21.1 days) was discovered serendip- 
itously by Liu et al. (|1992p in the course of spectro- 
scopic monitoring of SWCMa, which is only about 2. '5 
away. Beyond the initial estimates of the minimum masses 
by these authors (1.74Mq and 1.80 Mq), no determina- 
tions of the absolute dimensions have been reported in 
the literature. Accurate differential uvby photometry of 
HWCMa was published by Clausen et al. (|2008p . who 
confirmed that the combination of high orbital eccentric- 
ity (e « 0.50), low inclination {i « 84?7), and a line of 
apsides practically aligned with the line of sight result in 
the complete absence of eclipses near apastron. The shal- 
low (0.13 mag) eclipses that do take place near periastron 
correspond to the less massive star being occulted by the 
other. In the following we refer to these as the primary 
eclipses, following the convention adopted by Clausen et 
al. (|2008p . Secondary eclipses would be expected at phase 
~0.52, but as mentioned above, they do not occur. 

An accurate ephemeris has been determined by 
Clausen et al. (^005]) as 

Min I (HJD) = 2 452 279.6787(4) + 21^1178329(33) x E , 

which we adopt here. Apsidal motion is expected in view 
of the large eccentricity, but due to the lack of secondary 
eclipses this can only be measured spectroscopically. 

3.1. Radial velocities and spectroscopic orbit 

Spectroscopic observations of HW CMa were gathered be- 
tween 1988 November and 1999 March with the same 
instrumentation used for SWCMa and described in 
Sect. 12. ll A total of 48 spectra were included in the analy- 
sis, with S/N ratios ranging from 9 to 39 per resolution ele- 
ment of 8.5 kms~^. This is a superset of the observations 
reported by Liu et al. (|1992p . who used the same spec- 
trographs. The optimization of the synthetic templates 



and determination of radial velocities using TODCOR fol- 
lows the procedures described earlier. The surface grav- 
ities for the templates were held fixed at logg = 4.0 
for both stars, which is the value in our grid of syn- 
thetic spectra closest to the final estimates reported in 
Sect. 13.41 As was the case for SWCMa, cross-correlation 
grids with TODCOR indicated a preference for a metallic- 
ity of [Fe/H] = +0.50 for both stars, the highest available 
in our library of templates, but in this case we saw no ob- 
vious difference between the two stars. A value this high 
is again a hint of the possible Am nature of the com- 
ponents. The template parameters adopted for the RV 
determinations are Tcff = 7500 K (less massive star, the 
'primary') and Tog = 7750 K (more massive 'secondary'), 
and Dsini = 16 kms~^ for both components. 

Simulations with synthetic binary spectra were run in 
the same way as for SW CMa to estimate and correct for 
systematic errors in the raw velocities resulting from lines 
shifting in and out of our narrow spectral window with or- 
bital phase. The magnitude of these effects for HW CMa 
is less than 1 kms~^ (see Fig. [5]). Nevertheless, the im- 
pact on the masses and radii is not negligible in this case: 
the masses increase by about 1% (^1.5 times their final 
errors), and the radii by slightly more than 0.3%. The 
final velocities including these corrections are listed in 
Table lX^ and incorporate also the adjustments for instru- 
mental shifts from run to run described in Sect. 12.11 The 
RV measurements and our best-fit spectroscopic orbital 
solution are shown graphically in Fig. |9l and the elements 
are given in Table O The light ratio inferred from the 
CfA spectra using TODCOR is Lp/L^ = 0.92 ± 0.04 (pri- 
mary/secondary), which corresponds to the mean wave- 
length of those observations. 

The final surface gravities of the stars (logg « 4.24; 
Sect. 13. 4p are approximately halfway between two steps 
in our library of templates, and our velocity determina- 
tions above used the lower of the two values (logg = 4.0). 
Given that effective temperatures are strongly correlated 
with surface gravity in determinations based on our CfA 
spectra, a more accurate temperature may be derived 
by interpolation to logy = 4.24. In this way we find 
Teff = 7630 ± 150 K and 7780 ± 150 K for the primary and 
secondary, and v sin i values (including other sources of 
line broadening; see Sect. 13.41) of 15 kms~^ for both stars. 
We expect the temperatures to be more accurate than in 
the case of SW CMa, as the line strengths of the two stars 
in HW CMa appear more or less equally enhanced com- 
pared to solar metallicity stars. 

3.2. Chemical abundances 

For a detailed study of the chemical composition of 
HWCMa a spectrum of the binary was obtained on UT 
2008 February 26 with the FIES instrument on the 2.5 m 
Nordic Optical Telescope (La Palma). The Julian date at 
mid exposure is HJD 2 454 523.38203, corresponding to 
orbital phase 0.2469. The S/N ratio of this observation is 
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Table 5. Spectroscopic orbital solution for HWCMa. 



Table 6. Abundances for the components of HW CMa. 



Parameter 



Value 



Primary 



Secondary 



Adjusted quantities: 
Kp (kms-i) 
Ks (kms-^) 
7 (kms"''") 

Adopted quantities: 
P (days) 

Ti (HJD-2 400 000) ^ 
Derived quantities: 

Mp sin^j (Mq) 

Mssin^i (Mq) 

q = Mp/Ms 

ap sini (10^ km) 

as sin i [lO'' km) 

asini (Rq) 
Other quantities pertaining 

Nobs 

Time span (days) 
(jp (kms"'^) 
as (kms""*") 



68.47 ± 0.23 
66.18 ±0.21 
+16.91 ±0.10 
0.5016 ± 0.0018 
86.38 ± 0.26 

21.1178329 
52279.6787 

1.700 ±0.011 
1.759 ±0.012 
0.9665 ± 0.0043 
17.202 ± 0.055 
16.626 ± 0.050 
48.63 ±0.10 

to the fit: 

48 
3759 
1.05 
0.95 



^ Longitude of periastron for the less massive (primary) 
star. 

^ Time of central primary eclipse (eclipse of the less massive 
star). 



approximately 120 per pixel at a wavelength of 6160 A. 
The reduction and analysis were carried out in the same 
way as described in Sect. 12.21 For the VWA analysis we se- 
lected only lines with measured equivalent widths smaller 
than 65 mA in the composite spectrum (corresponding to 
about 140 mA for the primary and 120 mA for the sec- 
ondary), and we assumed a light ratio of Lp/Lg — 0.90 
for the V band (see Sect. 13. 4p . The surface gravities of 
the components were both held fixed at a preliminary es- 
timate of \ogg = 4.23, near the final values in Sect. 13.41 
The wsini values were measured to be 12.0 ± l.Okms"^ 
for both stars (for an adopted macroturbulent velocity of 
Crt = 7kms~^), although this has little impact on the 
abundance results. The temperatures yielding the best 
agreement between the abundances of Fei and Fell are 
7500 K for the primary and 7700 K for the secondary, with 
uncertainties estimated at 150 K. 

We tested a range of [Fe/H] values for the model at- 
mospheres used in the VWA analysis, as well as differ- 
ent values of the microturbulence in order to examine the 
sensitivity of the results. The final [Fe/H] level of +0.30 
dex for the model atmospheres was chosen to be consis- 
tent with that found below from the equivalent widths. 
Microturbulent velocities were 2.55 ±0.80 km s~^ for both 
stars. Individual elemental abundances for 16 species are 
presented in Table IHl on the scale of the solar abundances 



Element 


[X/H] 


a 


iV 


[X/H] 


a 


iV 


Cl 


-0.66 




1 








Nai 


+0.34 




1 


+0.37 




1 


Mgi 








-0.24 




2 


Si I 


+0.02 


0.09 


6 


+0.13 


0.12 


4 


Sin 


+0.40 




1 


+0.24 




2 


Cai 


-0.64 


0.17 


4 


-0.65 




2 


Tin 


+0.25 




2 


+0.18 


0.15 


4 


Or I 


+0.08 




1 


-0.31 




1 


Cm 


+0.57 


0.18 


3 


+0.35 


0.11 


6 


Fei 


+0.30 


0.08 


31 


+0.29 


0.08 


24 


Fell 


+0.37 


0.09 


12 


+0.35 


0.10 


15 


Nil 


+0.32 


0.15 


5 


+0.52 


0.09 


12 


Znl 


+0.28 




1 


+0.26 




2 


Sril 


+0.72 




1 


+0.10 




1 


Yii 


+0.87 




1 


+0.75 




1 


Ban 


+1.14 




2 


+1.16 


0.18 


3 



Note: These abundances are based on Vmic ~ 2.55 kms~^ for 
both stars, and [Fe/H] = +0.30 for the model atmospheres 
applied in the analysis. Fei abundances include corrections for 
NLTE. A'' represents the number of lines of each element, and 
a is the scatter of the measurements (listed only if A > 2). 
The primary is the less massive star in the binary. 



of Grevesse et al. (|2007p . An NLTE correction to Fei of 
+0.06 dex for the primary and +0.07 dex for the sec- 
ondary (Rentzsch-Holm I1996P is included. The average 
iron abundance is essentially identical for the two stars: 
[Fe/H] = +0.33 ±0.15 for the primary and +0.32 ±0.15 for 
the secondary (uncertainties include contributions from 
errors in Toff and Wmic)- 

The abundance pattern for the components of 
HW CMa is shown graphically in Fig. [TUl There are strong 
indications of the Am phenomenon in both stars: Ca seems 
characteristically deficient, along with C (at least in the 
primary), and heavier elements appear to be enhanced. 
In particular, Sr, Y, and Ba are generally much stronger 
than normal. The Fe/Ca difference is the same for both 
components (0.97 dex), and is considerably larger than 
we found for SWCMa. Additionally, for that system we 
measured Zn to be stronger than iron by a factor of 3-4 
in both stars, while in HW CMa these two elements have 
about the same abundance. Scandium is a highly diag- 
nostic element in Am stars, where it is usually found to 
be deficient. Unfortunately, however, no lines of this ele- 
ment are suitable for measurement in our spectra of ei- 
ther system. Nevertheless, the overall abundance pattern 
in both SWCMa and HWCMa is quite similar to that 
found in many other studies of Am stars (see, e.g., Okyudo 
& Sadakane IT9901 Hui-Bon-Hoa [20001 Yushchenko et al. 
Iliev et al. Fossati et al. [MTTl Gebran et al. 
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3.3. Photometric elements 

The uvby differential pliotometry used for HWCMa is 
tfiat reported by Clausen et al. (|2008p . obtained with 
the 0.5 m Stromgren Automatic Telescope at ESO (La 
Silla) between 1989 February and 2002 March. A total of 
415 observations were obtained in each filter, with similar 
average uncertainties per differential measurement as for 
SWCMa (~5 mmag in vby and ~6 mmag in u). The sys- 
tem was observed intensively at the expected phase of sec- 
ondary eclipse, but no significant drop in brightness was 
detected, as mentioned earlier. With only a single eclipse 
to constrain the geometry, this data set therefore poses 
special challenges to extract reliable orbital elements: the 
central surface brightness ratio Js is completely uncon- 
strained, there is no information in the light curve on the 
eccentricity (ecosa;), and the radius ratio k is very poorly 
defined without an external constraint. 

As for the case of SW CMa, we used the JKTEBOP 
program for the light-curve solutions given that the stars 
are essentially spherical (oblateness ~ 0.00006). The 
ephemeris adopted is that of Sect. [31 The eccentricity 
and longitude of periastron were held fixed at the spec- 
troscopic values, as was the mass ratio. In order to con- 
strain Js we calculated model spectra for each star from 
the MARCS series (Gustafsson et al. '20081 for tl^e adopted 
temperatures (7560 K and 7700 K for the primary and sec- 
ondary; see Sect. 13. 4p . and used them to obtain Js esti- 
mates of 1.078 (u), 1.105 (v), 1.091 (6), 1.075 (y), and 
1.081 (5187 A), with only a small dependence on [Fe/H] 
and logg. The y value is nearly identical to that obtained 
from the empirical flux scale of Popper (|1980[) . Js = 1.076 
in V, supporting the accuracy of the theoretical calcula- 
tions. 

Because of its sensitivity to k, the spectroscopic light 
ratio [Ls/Lp oc Jsk^) provides a very useful constraint 
on the poorly defined ratio of the radii in HWCMa. The 
average spectroscopic light ratio from three independent 
determinations described below in Sect. 13.41 is Lp/ Ls = 
0.90 ± 0.02, or Ls/Lp = 1.11 ± 0.02. This corresponds 
strictly to the mean wavelength of our spectroscopic ob- 
servations (5187 A). Light ratios in the Stromgren bands 
that are needed for the light curve fits were inferred by 
making an initial estimate of k requiring that the combi- 
nation of the measured Lg/ Lp and k match the theoreti- 
cal value of Js at 5187 A. We then used this estimate of 
k = 1.013 along with the Js values in uvby to obtain the 
corresponding hght ratios: 1.106 (u), 1.134 (d), 1.119 (6), 
and 1.103 (y). 

JKTEBOP fits were performed separately in each pass- 
band, fixing Js and Ls/Lp to the values specified above 
and holding e sin u and e cos w fixed as well (from the spec- 
troscopy). We solved for i, Vp + Ts, k, and the usual pho- 
tometric scale factor and phase offset. Solutions were car- 
ried out using LD coefhcients for the linear law from both 
the Claret (l2000l) and van Hamme (fT993)) tables. The ge- 
ometric elements with each set of coefficients show good 
agreement between the v, 5, and y bands, with u being 



Table 7. Adopted photometric elements for HW CMa. 



in 


84.84 ± 0.06 


rp + rs 


0.06769 ± 0.00066 


k 


1.012 ±0.013 




0.03365 ± 0.00037 


rs 


0.03404 ± 0.00042 





y 


b 


V 


u 


Js 


1.075 


1.091 


1.105 


1.078 




±2 


±2 


±2 


±2 


Ls / Lp 


1.100 


1.121 


1.135 


1.108 


±20 


±20 


±21 


±20 



Note: Flux ratios are based on the measured eflfective tem- 
peratures and MARCS models. Luminosity ratios are con- 
strained to match the spectroscopic value at 5187 A. The geo- 
metric elements are the weighted mean of the vby values from 
Table IC.2I using hnear limb-darkening coefficients from van 
Hamme (|T993]l . 



more discrepant (and also more uncertain). We describe 
the results of these tests in Appendix [C] (available elec- 
tronically) . 

In previous papers of this series the LD coefficients 
from van Hamme (1993) have generally been found to be 
in better agreement than those from Claret (|2000p with 
the values that result when the coefficients are left free 
in the solutions (which cannot be done here, for obvious 
reasons), and that was the case for SWCMa as well (see 
Appendix |B|) . Consequently, for the final light elements 
of HW CMa we adopt the van Hamme (1993) coefficients, 
and take the weighted average of the vby solutions (ex- 
cluding u). These values are presented in the top portion 
of Table [71 where the uncertainties include a contribution 
from the scatter between the different bandpasses. A final 
set of solutions was carried out in each passband with the 
geometric elements held fixed at these weighted average 
values, in order to derive the luminosity ratios. These are 
listed at the bottom of the table. The observations and fit- 
ted curves near the primary eclipse are shown in Fig. 1111 
The final fits indicate only 22% of the y-band light of 
the photometric primary (less massive star) is blocked at 
phase 0.0. 

3.4. Absolute dimensions 

The mass determinations for HW CMa are both good to 
about 0.7%. Despite the difficulties described in the pre- 
ceding section, the radii are also very precise, with relative 
uncertainties of only 1.1% and 1.3% for the primary (less 
massive star) and secondary that we consider realistic. 
These and other properties of the system are collected in 
Table m 

With the uvby/3 indices out of eclipse from Clausen 
et al. (|2008p and the calibration by Crawford (|1979p we 
derive an estimate of the interstellar reddening toward 
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HW CMa of E{b - y) = 0.026. We assign to this a conser- 
vative uncertainty of 0.01 mag. The corresponding extinc- 
tion is A{V) = 0.11 ± 0.04. A rather similar value is ob- 
tained from the dust maps of Hakkila et al. (|1997p , which 
give A{V) — 0.13 or E{b — y) = 0.031 after correction for 
the distance to the binary. On the other hand, Schlegel 
et al. (|1998p suggest larger values of A{V) = 0.20 and 
E{b — y)= 0.049. In the following we adopt the more reli- 
able photometric estimate. Table [5] includes also the indi- 
vidual Stromgren indices of the components, derived from 
the combined uvby magnitudes and the light ratios used 
in our photometric solutions. After de-reddening, these 
lead to photometric temperatures of 7777 ± 102 K and 
7972 ± 107 K using the calibration of Napiwotzki et al. 
(|1993p . in which the uncertainties include observational 
errors but not the unknown dispersion of the calibration. 



Table 8. Astrophysical properties of the HWCMa sys- 
tem. 





Primary 


Secondary 


Absolute dimensions: 






M {Mq) 


1.721 ±0.011 


1.781 ±0.012 




1.643 ±0.018 


1.662 ±0.021 


log 3 (cgs) 


4.242 ±0.010 


4.247 ±0.011 


V sin j (km s~^) 


12.0 ± 1.0 


12.0 ± 1.0 


Vsync (kms^^) 


3.9 ±0.1 


4.0 ±0.1 


«psync (kmS"^) " 


11.1 ± 0.1 


11.2 ±0.2 


Dperi (kmS"^) 


13.7 ±0.2 


13.8 ±0.2 


Radiative and other properties: 




TeS (K) 


7560 ± 150 


7700 ± 150 


Mbol (mag) 


2.488 ± 0.089 


2.383 ±0.088 


log L/Lq 


0.898 ± 0.036 


0.940 ±0.035 


BCv (mag) 


±0.03 ±0.10 


±0.03 ±0.10 


Mv (mag) 


2.46 ±0.13 


2.35 ±0.13 


Distance (pc) 


306 ± 15 


V - Mv (mag) 


7.43 ±0.10 


[Fe/H] <= 


±0.33 ±0.10 


±0.28 ±0.10 


Photometric indices: 








9.996 ±0.012 


9.892 ±0.011 




0.139 ±0.005 


0.118 ±0.005 


mi * 


0.225 ±0.011 


0.231 ±0.011 




0.817 ±0.016 


0.858 ±0.016 


E{b-y) 


0.026 ±0.010 



^ Measured projected rotational velocity. 
^ V sin i expected for synchronous rotation. 

vsini expected for pseudo-synchronous rotation. 

Expected v sin i if synchronized at periastron. 

Values representative of the surface layers only. 
' Not corrected for interstellar absorption/reddening. 

Note: The bolometric corrections (BCv) are adopted from 
Flower (1996), along with with T^^q = 5780 K, and M^^iq = 
4.73 (see Torres I2010[l . An additional error contribution of 
0.10 mag is added in quadrature to the BCv uncertainty prop- 
agated from the temperature uncertainties. 



These estimates are 200-250 K hotter than the adopted 
spectroscopic values. 

In addition to the two consistent estimates of the effec- 
tive temperatures reported earlier from the CfA/DS and 
FIES spectra, we obtained another measure by disentan- 
gling the DS spectra and using a cross-correlation proce- 
dure similar to that described in Sect. 12.11 for SWCMa, 
but adapted for single-lined spectra. The disentangling 
was carried out with a revised version provided by E. 
Sturm of the original procedure introduced by Simon 
& Sturm P^M)) . The resuhs are 7500 K and 7620 K for 
the primary and secondary^ A further spectroscopic es- 
timate of the temperatures was made from three spec- 
tra of HWCMa collected with the TRES instrument 
(Szentgyorgyi & Friresz [2007p on the 1.5 m reflector at 
the F. L. Whipple Observatory (A/AA « 44 000, AA3900- 
8900 A). These observations were initially intended to 
support the abundance determinations in Sect. 13.21 but 
turned out to be too weak for that purpose. Application 
of the TODCOR-based cross-correlation procedures used 
before for composite spectra gave Toff values of 7610 K 
and 7670 K for the primary and secondary, and a light ra- 
tio oiLp/Ls = 0.89±0.02 at a mean wavelength of 5187 A. 
Finally, we applied the same TODCOR analysis technique 
to the FIES spectrum we used earlier for the abundance 
analysis, and obtained 7570 K and 7740 K along with an 
identical light ratio of 0.89 ± 0.02. The above determina- 
tions are all in good agreement, and establish that the pri- 
mary (less massive star) is slightly but significantly cooler 
than the secondary. We adopt a straight average of the 
five temperature estimates, 7560 K and 7700 K for the pri- 
mary and secondary, to which we attach an uncertainty of 
150 K. The corresponding spectral type is approximately 
A6. The average of the three independent spectroscopic 
light ratio determinations is Lp/Lg = 0.90 ± 0.02. 

The luminosities and distance were derived in the same 
way as for SWCMa. The system distance is 306 ± 15 pc. 
Nearly identical values of 307 ± 20 pc and 306 ± 20 pc ob- 
tained separately for the primary and secondary indicate 
excellent internal consistency. HW CMa and SW CMa are 
only 2. '5 apart on the sky and HW CMa is only half as dis- 
tant as SW CMa, yet the reddening we find for HW CMa 
is slightly larger. This may indicate that interstellar ex- 
tinction is patchy in this direction. 

The measured projected rotational velocities for both 
components of HW CMa are consistent with the pseudo- 
synchronous velocities, assuming co-aligned orbital and 
spin axes (see Table [5]). 



^ While the cross-correlation technique applied to these dis- 
entangled spectra is able to extract a useful estimate of Tes 
because it makes use of the full spectrum at once, there are not 
enough isolated lines of iron and other elements in this narrow 
wavelength region to attempt a detailed abundance analysis as 
was done earlier with the FIES spectrum. 
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4. Discussion 

The following sections present a comparison of the mea- 
sured properties for SWCMa and HWCMa against pre- 
dictions from the theory of stellar evolution, stellar struc- 
ture, and tidal evolution. 

4.1. Stellar evolution 

Accurate mass, radius, and effective temperature determi- 
nations, especially when accompanied by a measurement 
of the chemical composition, are among the most powerful 
constraints afforded by eclipsing binary stars to test mod- 
els of stellar structure and evolution (see, e.g., Andersen 
119911 Torres et al. 120101 and references therein). When 
these four quantities are available for both components 
in the system, there are essentially no free parameters in 
the comparison with publicly available models: to the ex- 
tent that the measurements are accurate, the models will 
either be successful in reproducing all observations simul- 
taneously within their errors, or they will failH 

In the case of Am stars such as SW CMa and HW CMa, 
however, the measured abundances do not represent the 
bulk composition of the stars, which is the relevant quan- 
tity when comparing with models. Instead, they are only a 
reflection of local changes in the surface layers attributed 
to diffusion processes, which give the spectra of these ob- 
jects their peculiar appearance. Metallicity therefore re- 
mains a free parameter when testing models for these 
stars. Nevertheless, because detailed abundance patterns 
in Am stars vary from case to case and are not yet com- 
pletely understood, a comparison with models is still of 
considerable interest ~ even if less constraining ~ to es- 
tablish the precise evolutionary state of the stars. Except 
for second-order effects due to limb darkening, absolute 
masses and radii are essentially independent of the the- 
ory being tested. For HW CMa, however, our light-curve 
solutions have had to make additional use of model atmo- 
spheres to constrain the geometry of the system, because 
of the lack of secondary eclipses. We point out, though, 
that those models were used only in a differential sense 
(to determine ratios of central surface fluxes or bandpass- 
specific luminosities), so in effect the dependence on the- 
ory is very weak. 

In this section we compare the absolute dimensions 
of SW CMa and HW CMa with stellar evolution models 
from the Yonsei-Yale series (Yi et al. 120011 Demarque et 
al. I2004P , which have been shown previously to match ob- 
servations very well. We have chosen these models primar- 
ily because they are provided with software to interpolate 
evolutionary tracks (and isochroncs) for any given mass 

® Additional variables such as the helium abundance, mix- 
ing length parameter, convective core overshooting, etc., have 
generally already been chosen in advance by the modelers in 
published tables, and are not "tunable" by the user. Specific 
models for a particular system can of course be computed with 
different values of these additional variables, as we do in one 
example below, adding more freedom to the fits. 



and metallicity, the mass being among the most accurately 
known properties for these stars. 

The log g vs. Toff diagram for SW CMa in Fig. [H al- 
lows a simultaneous comparison of the models against 
the three key observables: M, R, and Teff. Evolutionary 
tracks corresponding to a composition slightly above so- 
lar ([Fe/H] = -1-0.05, with no a-element enhancement) 
are shown for the exact masses we measure (Table |4]), 
with the gray areas representing the uncertainty in the 
location of the tracks that comes from the mass errors 
(in this small fraction of the mass difference). 

Isochroncs are also displayed for ages 500-900 Myr, in 
steps of 100 Myr. The models for this composition pro- 
vide a very good fit to the observations, within the uncer- 
tainties, and the age inferred is approximately 700 Myr 
(indicated with the thicker isochrone). SW CMa is seen to 
be somewhat evolved, and currently near the mid-point of 
its main sequence phase. The implied bulk composition in- 
ferred from theory, [Fe/H] = -1-0.05, is ~0.4-0.5 dex lower 
than the iron abundance measured at the surface ([Fe/H] 
= -I-0.49/-I-0.61 for the primary and secondary, respec- 
tively), which may be interpreted as giving a rough mea- 
sure of the true enhancement of the iron-peak elements 
in the photospheres of these stars. Strictly speaking, how- 
ever, this factor of three enhancement may only be a lower 
limit, as the effects of metallicity and a-element enhance- 
ment in the models can trade off against each other to 
some extent, leading to mass tracks and isochroncs that 
fit the observations equally well. For example, a match as 
good as seen above is achieved also with [Fe/H] = —0.10 
and [a/Fe] = +0.22, with [Fe/H] = -0.20 and [a/Fe] 
= +0.36, etc. 

Separate diagrams of the radius and effective temper- 
ature as a function of mass are seen in Fig. [131 along with 
the same isochroncs from the previous figure. The models 
provide a virtually exact match to the mass and radius of 
the components at 700 Myr, and a good fit to the temper- 
atures as well. 

Similar diagrams for HW CMa can be seen in Fig. [T3] 
and Fig. [151 In this case, the model composition that 
best fits the observations is [Fe/H] = +0.23 (with [a/Fe] 
= 0.0), and the implied age is approximately 160 Myr. 
Solar composition tracks are shown for reference as dot- 
ted lines in Fig. [TH HW CMa is essentially on the zero-age 
main sequence. To the extent that [Fe/H] = +0.23 rep- 
resents the overall composition of this system, it would 
appear that the true enhancement of the iron-peak ele- 
ments in the surface layers is a modest 0.1 dex in this case, 
since the measured photospheric abundances are [Fe/H] 
= +0.33/+0.32. As before, this conclusion is dependent 
on the assumed [a/Fe] enhancement, which is inaccessi- 
ble to observation in both of these systems. The higher- 
than-solar bulk [Fe/H] = +0.23 composition for HWCMa 
suggested by the models is somewhat uncommon for stars 
in the solar neighborhood. If one enforces in the models 
a more typical solar-like iron abundance, the observations 
are about as well reproduced by increasing the [a/Fe] ratio 
to +0.40. This combination would imply a true surface- 
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layer iron enhancement due to the Am phenomenon in 
HW CMa of '--^0.3 dex, and the inferred age would be the 
same as before. However, a value of [a/Fe] as large as 
this is on the high side of the observed distribution for 
stars with solar iron abundance, which makes us some- 
what skeptical of this scenario. Another possible way of 
reproducing the HW CMa observations with an iron com- 
position closer to solar is to alter the helium abundance. 
Tests with the Granada models of Claret (j2004p indicate, 
however, that the radii and effective temperatures can 
only be matched simultaneously with a very low helium 
abundance of 1" = 0.24, at ages of 440-500 Myr that are 
significantly older than before. A helium abundance this 
close to the primordial value seems unlikely. 

Comparisons (not shown) of the observations for 
HWCMa and SWCMa with the standard ([a/Fe] = 0.0) 
Victoria- Regina stellar evolution models of VandenBerg et 
al. (|2006D and with the Granada models of Claret (|2004p 
give results that are quite consistent with those from the 
Yonsei-Yale models, indicating similar overall abundances 
and evolutionary ages. 

4.2. Internal structure 

The observed apsidal motion for SWCMa, Co = 0.00067 ± 
0.00021 deg cycle^^ (Clausen et al. I2008p . is a measure 
of the average degree of central mass concentration of 
the components. The corresponding average apsidal mo- 
tion constant for the system is logA:2,obs — —2.57 ± 0.30. 
Using the absolute dimensions of the stars from Table 21 
the models by Claret (|2004p predict a theoretical value 
of logfc2.thco = —2.582 ± 0.050, in which the formal un- 
certainty does not account for possible systematic errors 
in the models. While this is consistent with the empirical 
value within the errors, the close agreement may be acci- 
dental given the difficulty of the measurement (the apsi- 
dal period is long: U = 14 900 ± 4700 yr). The General 
relativistic contribution (wgr = 0.00034 deg cycle^^) 
amounts to about 50% of the total apsidal motion. 

While the eccentric system HW CMa is also expected 
to display apsidal motion, the lack of secondary minima 
precludes such a measurement based on eclipse timings. 
Spectroscopic detection of the effect would be quite chal- 
lenging. The models of Claret (|2004p predict a value of 
uj = 0.00023 deg cycle"^, and a corresponding apsidal 
motion constant of logfc2,thco — —2.410 ± 0.050. The es- 
timated period of this effect is U ~ 88 900yr, six times 
longer than in SW CMa. In the case of HW CMa the rela- 
tivistic contribution completely dominates, accounting for 
97% of the total apsidal motion. 

4.3. Tidal evolution 

In addition to their binarity. Am stars typically display 
slow rotation compared to field stars of the same spec- 
tral type. The components of SW CMa and HW CMa are 
no exception. We noted earlier that while the primary of 



SW CMa appears to be pseudo-synchronized with the or- 
bital motion, the secondary is not, and rotates at only 
half the corresponding speed. This assumes the spin axes 
are parallel to the axis of the orbit, a condition that is 
almost universally taken for granted in binaries. Here we 
investigate whether synchronization is expected at all in 
these systems according to tidal theory, and if so, whether 
the situation for the secondary of SW CMa might be un- 
derstood in terms of spin-orbit misalignment. We examine 
also the predicted evolution of the orbital eccentricity. 

We have integrated the differential equations given 
by Hut (|198ip for the time evolution of the semimajor 
axis {da/dt), the eccentricity {de/dt), the angular rota- 
tion rates of both components {dfti/dt, dil2/dt), and the 
angle between the plane of the orbit and the equator of 
each star {d(j)i/dt, d(j)2/dt). These six, coupled equations 
were integrated simultaneously using a 4th-order Runge- 
Kutta method, with the stellar properties interpolated at 
each time step from appropriate evolutionary tracks from 
Claret (|2004p . The turbulent dissipation timescale for the 
late evolutionary phases with convective envelopes was 
taken to be (Mi?^/^)^/^, where M, R, and L are the 
mass, radius, and luminosity of the star. For phases in 
which the envelopes are radiative (the most relevant here) 
the timescales adopted follow closely those in Eq. 17 and 
Eq. 18 by Claret & Cunha (fTWT)) . 

The initial conditions are of course not known, so they 
are in effect free parameters of the model. For the orbital 
period and eccentricity we have adjusted the starting val- 
ues so as to match P and e at the current evolutionary age 
of each system. Their time evolution for SW CMa is shown 
in Fig.[T5k and Fig.lTBb. in which the vertical dotted lines 
mark the current age of about 680 Myr (logr = 8.833), 
according to the Granada models employed in this sec- 
tion. Based on these calculations the orbit is expected to 
circularize at an age of r = 835 Myr (logr — 8.921) for 
this system. 

The evolution of the angular rotation rate f2 = 27r/Pi-ot 
of each star in SWCMa is seen in Fig. [TBh . where for 
convenience we have normalized it to the orbital rate, 
r^orb = 27r/Porb- Wc havc no observational constraint 
on the rotation rates, so the initial values were arbitrar- 
ily set to r2i/J7orb = f^2/f^orb = 6.0, reasonable for A- 
type stars that usually rotate very rapidly (at least in 
the field). A dot-dashed line in this figure indicates the 
pseudo-synchronous rate at each age, and the horizontal 
dotted line is the current pseudo-synchronous rate. These 
calculations suggest that pseudo-synchronization for both 
stars in this system was reached early-on, at an age of 
roughly 100 Myr (logr = 8.0). This would indeed appear 
to be the case for the primary of SW CMa, based on its 
measured vsini value, but not for the secondary. 

The degree of spin-orbit alignment is represented by 
the angle 0, but once again, we have no constraint on 
this quantity. Furthermore, the relation between cj) and 
the orbital and rotational inclination angles iorb and irot, 
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at the current age of the system, (f) Same as (e), for the secondary. In this case theory is not able to match the 
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both measured with respect to the line of sight, is given 

by 



cos 4> = cos iorb COS irot + sin Zorb sin irot COS A 



(1) 



which involves an unknown angle A between the sky- 
projected angular momentum vectors of the orbit and the 
stellar spin. 

The spectroscopically measured projected rotational 
velocities of the stars, which we refer to more properly 
now as V sin Zj-ot , do provide an indirect constraint on a 
combination of theoretically predictable quantities, but 
this still involves the unknown angle A. Given that iorb 
is rather close to 90° in both SW CMa and HW CMa we 
may make the approximation that cos (p Ri sin ij-ot cos A. 
We then have 



27r rt cos (j) 
tismirot ~ "5 — r-R • 

^orb ^ ^orb COS A 



(2) 



All quantities on the right-hand side of this equation are 
either known from stellar evolution calculations (i?), or 
can be computed from the solution of the differential 
equations for tidal evolution, with the exception of the 
angle A, which depends on the observer's viewpoint. In 
order to make progress, we ignore this term for the mo- 



ment (or equivalently, we consider A to be small), so that 
cos (f) w sin ii-ot . 

The evolution of the alignment angle cf> for SW CMa is 
displayed in Fig. [T5H for four different initial values (20°, 
40°, 60°, 80°). The curves for the primary and secondary 
are nearly indistinguishable, so only those for the primary 
are shown. As seen from the convergence of the curves to- 
wards zero, spin-orbit alignment for both stars in SW CMa 
is reached at the age of about 200 Myr (logr « 8.3), which 
is less than a third of the current age of the system. The 
predicted evolution of wsinij-ot for each star, computed 
from Eq.[51 is shown in Fig. \Wk and Fig. [TBF compared 
against the measured values (dashed line and shaded un- 
certainty region). 

Independently of the initial values of (j) and fi/fiorb, the 
predictions from tidal evolution theory cannot be made 
to match the anomalously slow rotation rate of the sec- 
ondary of SW CMa. The angle A in Eq.[2] that we have 
previously ignored is inconsequential for this discussion, 
as any value greater than zero only makes the disagree- 
ment worse. Evidently current tidal theory is incomplete. 
An effect not considered here is the possible decoupling be- 
tween the core of the star and the external layers, which 
suffer more directly the actions of tidal forces. Since the 
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spectroscopically measured v sin ij-ot corresponds only to 
the layers that we see, it is quite possible that the larger 
rotation rate predicted by theory is more representative of 
a more rapidly rotating core than of the outer envelope, or 
at the very least is some average of the two, which would 
be larger than the measured value. 

Similar tidal calculations for HWCMa are shown in 
Fig-HZl In this case we have chosen to illustrate the effects 
of a change in the initial values of r2/f2orb, rather than (f>. 
As before, we have adjusted the boundary conditions for 
the period and eccentricity to match the measured val- 
ues at the current age, which is 160 Myr (logr — 8.21) 
according to the Claret (|2004p models. As the evolution 
of both of these quantities (Fig. [T7k and Fig. \T7h) de- 
pends on Q (see Hut 11981]) . the starting values for P and 
e are different in each case. The calculations indicate that 
orbital circularization will occur at an age of ~1.75Gyr 
(logr = 9.244). 

The time dependence of 0/17orb appears in Fig. [T7b . 
and influences the late evolution of the angle (j) between 
the equator and orbital plane, as seen in Fig. [T7b . These 
two panels indicate that neither star is expected to be 
pseudo-synchronized at the current age, nor has spin- 
orbit alignment been achieved. Pseudo-synchronization is 
predicted to happen at ^600 Myr (logr — 8.78), and 
alignment at an age of approximately 1 Gyr. The fact 
that we do see agreement within the measurement errors 
between the observed rotational velocities and the pre- 
dicted pseudo-synchronous velocities at the current epoch 
as listed in Table [5] suggests that this may be due to a 



combination of non-synchronized angular rotation rates 
{^1,2) and slightly misaligned spin axes ((^1,2 ~ 20° at 
the current epoch for our particular choice of the initial 
values for ^1^2; see Fig. [TTli). The theoretical projected 
rotational velocities are plotted in Fig. [T7b and Fig. [T7t . 
The lower curves match the measured v sin ij-ot of each star 
at the current age, and correspond to initial fJ/rJorb values 
of 4.0 for both stars, and initial angles 0i_2 — 30° tuned 
to produce the agreement. 

The asynchronous rotation of the secondary of 
SWCMa is not unusual for Am stars, or for normal A 
stars for that matter. The compilation of Torres et al. 
(2010) includes several such systems with well determined 
properties that show the same lack of synchronization. 
An extreme case is V459 Cas (Lacy I2004p . a pair of 
Aim stars both rotating 3.5-4.5 times faster than pseudo- 
synchronous. There is no obvious correlation that we can 
see between asynchronous rotation and other global prop- 
erties of Am or normal A stars. This issue has also been 
studied by others (Budaj 1996] [TM71 Ihev et al. 1^0051 
Prieur et al. I2006P , but no credible pattern has emerged. 

5. Final remarks 

SW CMa and HW CMa join the very small ranks of Am 
stars with accurately known absolute dimensions (better 
than 2% relative errors in M and R) that also have their 
detailed abundances determined, in this case for a dozen 
or more elements in each star. We find that the measured 
masses, radii, and temperatures are well matched by stel- 
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lar evolution models such as those by Yi et al. (|2001l) . 
with inferred bulk abundances of [Fe/H] = +0.05 and 
+0.23 and ages of ~700 and '^IGOMyr, respectively. These 
abundances are 0.4-0.5 dex and 0.1 dex lower than those 
measured in the outer layers. These two systems seem to 
confirm previous evidence that the properties of Am stars 
are generally matched by stellar evolution models just as 
well as those of normal A stars, suggesting their prop- 
erties are not fundamentally different, even though their 
surface abundances are. There are exceptions to this rule, 
though, such as the recently studied Am system of XY Get 
(Southworth et al. l201ip . whose stars do not seem to agree 
well with theory. 

Internal structure models predict that apsidal motion 
in both SWCMa and HWCMa (empirically determined 
only in the first case) should be dominated by General 
relativistic effects 50% and 97% contributions, respec- 
tively). Current tidal theory suggests that the angular 
rotation rates of the SWCMa stars should be pseudo- 
synchronized with the orbital motion by now, and that 
their spin axes should be parallel to the orbital axis. 
Indeed we find that the measured projected rotational ve- 
locity of the primary agrees with expectations, but the 
measured v sin i of the secondary is too small by a factor 
of two, even though its mass is only 6% smaller than the 
primary. This disagreement cannot be resolved by tuning 
free parameters in the models, and underscores our in- 
complete understanding of these processes. For HW CMa 
theory predicts that neither pseudo-synchronization nor 
spin-orbit alignment have yet been reached. However, the 
measured rotations do agree with the pseudo-synchronous 
values in this case, which can nevertheless be explained 
with a proper combination of initial values for il/riorb and 
the spin-orbit alignment angle 0. 
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Fig. 1. Main sequence eclipsing binaries taken from Torres 
et al. (|2010p in which both components have masses in the 
range 1.6-2.4 Mq, and relative errors in the masses and 
radii of 3% or better. Primary and secondary components 
are connected with dotted lines, and stars that are con- 
firmed or probable Am stars are marked with open circles. 
SW CMa and HW CMa are shown with filled symbols. The 
thick solid lines are drawn for reference, and correspond 
to solar-metallicity isochrones from the series of Yi et al. 
(|2001|) for ages of 50 Myr and 1 Gyr. 
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Radial velocity (km s~') 

Fig. 2. Systematic errors to the raw TODCOR veloci- 
ties of SWCMa determined from simulations based on 

synthetic binary spectra. Filled circles correspond to the 
primary, and open circles to the secondary. The velocity 
differences are shown as a function of orbital phase and 
radial velocity, and have been applied to the raw veloci- 
ties as corrections. Phase 0.0 corresponds to the primary 
eclipse. 
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Fig. 3. Measured radial velocities for SWCMa and or- 
bital solution (solid line for the primary, dashed for the 
secondary). The dotted line in the top panel indicates the 
center-of-mass velocity of the binary, and phase 0.0 corre- 
sponds to the primary eclipse. Residuals are shown at the 
bottom. 
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Fig. 4. Elemental abundances measured for SW CMa. 
Dark gray areas correspond to the primary, and light gray 
to the secondary. Filled and open symbols represent neu- 
tral and ionized species, respectively. 
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Fig. 5. uvby photometry for SWCMa together with our Fig. 7. Enlarged view of Fig. [5] around the secondary min- 
best fit models (Table [3]). Phase 0.0 corresponds to imum of SWCMa. 
the eclipse of the more massive star (primary eclipse). 
Residuals are shown at the bottom. 
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Fig. 6. Enlarged view of Fig. [5] around the primary mini- 
mum of SW CMa. 
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Fig. 8. Systematic errors in tlie raw TODCOR velocities 
of HW CMa determined from simulations based on syn- 
thetic binary spectra. Filled circles correspond to the more 
massive star, and open circles to the other component. 
The velocity differences are shown as a function of orbital 
phase and radial velocity, and have been applied to the 
raw velocities as corrections. Phase 0.0 corresponds to the 
eclipse of the less massive star. 
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Fig. 9. Measured radial velocities for HW CMa and or- 
bital solution (solid line for the more massive 'secondary' 
star, dashed for the primary). The dotted line in the top 
panel indicates the center-of-mass velocity of the binary, 
and phase 0.0 corresponds to eclipse of the less massive 
component. Residuals are shown at the bottom. 




Fig. 10. Elemental abundances measured for HWCMa. 

Dark gray areas correspond to the primary (less mas- 
sive star), and light gray to the secondary. Filled and 
open symbols represent neutral and ionized species, re- 
spectively. 
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Fig. 11. uwfey photometry for HWCMa near the prhxiary 
minimum, together with our best fit models (Table [?])• 
Residuals are shown at the bottom. 
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Fig. 12. Measurements for SWCMa compared against 
Yonsei-Yale models by Yi et al. (PM]) for [Fe/H] = +0.05 
and [a/Fe] = 0.0. Evolutionary tracks for the measured 
masses are shown with solid lines and shaded areas indi- 
cating the uncertainty in the location of each track coming 
from the mass errors. Isochrones from 500 Myr to 900 Myr 
in steps of 100 Myr are represented by the dashed lines. 
The best fitting 700 Myr isochrone is drawn with a thicker 
line. 
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Fig. 13. Effective temperatures and radii for SWCMa 
shown as a function of the measured mass, with the obser- 
vational errors represented by the shaded boxes. They are 
compared against [Fe/H] = +0.05 Yonsei-Yale isochrones 
by Yi et al. (POOT]) from 500 Myr to 900 Myr, in steps 
of 100 Myr, with [a/Fe] = 0.0. The best- fitting model is 
indicated with a heavier line. 
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Fig. 14. Similar to Fig. [HI but for HW CMa. The best- 
fitting metallicity is [Fe/H] — +0.23. Solar metallicity 
tracks are shown for reference (dotted lines). The dashed 
lines represent Yonsei-Yale isochrones from 100 Myr to 
1 Gyr in steps of 100 Myr. The best fitting 160 Myr 
isochrone is drawn with a thick solid line. 
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Fig. 15. Similar to Fig.Hl but for HWCMa. In this case 
the best-fitting metalhcity is [Fe/H] = +0.23. The dotted 
lines represent Yonsei-Yale isochrones from 100 Myr to 
1 Gyr in steps of 100 Myr. The best fitting 160 Myr model 
is drawn with a heavier solid line. 
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Table A.l. Radial velocities for SWCMa (with all cor- 
rections applied) and residuals from the final spectroscopic 
orbit. 
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Table A. 2. Radial velocities for HWCMa (with all cor- 
rections applied) and residuals from the final spectroscopic 
orbit. 
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0.74 


0.82 


1.13 


0.27 


49316.8629 


0.7008 


-16.18 


50.58 


1.03 


1.14 


-0.85 


0.30 


49330.7929 


0.3604 


40.55 


-7.91 


1.12 


1.24 


-0.37 


0.02 


49400.6145 


0.6667 


-9.17 


44.14 


0.89 


0.99 


-0.02 


0.26 


50003.0181 


0.1925 


70.08 


-37.97 


0.48 


0.53 


0.25 


-0.14 


50027.8831 


0.3699 


37.97 


-6.92 


0.96 


1.06 


-1.39 


-0.61 


50030.9707 


0.5162 


17.50 


17.33 


0.86 


0.95 


1.54 


-0.57 


50033.9660 


0.6580 


-8.13 


41.78 


0.82 


0.91 


-0.52 


-0.51 


50034.9653 


0.7053 


-16.29 


51.75 


0.61 


0.68 


-0.12 


0.60 


50059.8947 


0.8858 


-50.61 


86.64 


0.55 


0.61 


-0.44 


0.32 


50088.8617 


0.2575 


58.02 


-25.57 


0.50 


0.55 


-0.38 


0.43 


50092.8003 


0.4440 


27.65 


5.77 


0.88 


0.98 


0.20 


-0.25 


51239.7317 


0.7550 


-24.38 


61.40 


0.87 


0.97 


1.33 


0.38 
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Appendix B: Light curve solutions for SWCMa 
under different limb-darkening assumptions 

Tables HIT] and MM present our JKTEBOP solutions for 
fixed linear LD coefficients (up, Us) from van Hamme 
(|1993p and from Claret (12000p . respectively. The uncer- 
tainties we report are formal errors from the iterative 
least-squares procedure, which are often found to be un- 
derestimated. In each case the agreement between the 
uvby passbands is quite good, and the results using the 
different LD coefficients show small but perhaps signifi- 
cant systematic differences, with the inclination angle be- 
ing about 0?22 smaller using the Claret values, the sum 
of the radii being ^0.5% larger, and k also being 1-3% 
larger (resulting in values up to 1% smaller, and val- 
ues 1-3% larger). Experiments using Claret coefficients for 
a quadratic law instead of the linear law lead to solutions 
that are not significantly better and results that are not 
very different, and tend to be closer to those obtained with 
the van Hamme coefficients. 

Solutions with the LD coefficients free are presented 
in Table lB73l subject only to the condition that the coef- 
ficients be the same for the primary and secondary since 
their temperatures are also very similar. These results in- 
dicate slightly better agreement than before between the 
four passbands. On average the fitted LD coefficients are 
closer to those by van Hamme than those by Claret. We 
summarize the elements obtained from the three differ- 

Table B.l. Photometric solutions for SWCMa with lin- 
ear limb darkening coefficients adopted from van Hamme 
(fT993)l . 



Parameter 


y 


b 


V 


u 


i n 


88.58 


88.71 


88.75 


88.53 


±5 


±3 


±3 


±7 


ecosctj 


-0.30399 


-0.30314 


-0.30318 


-0.30318 




±5 


±5 


±5 


±8 


esintj 


+0.0906 


±0.0880 


±0.0877 


±0.0900 




±13 


±12 


±12 


±20 




0.3163 


0.3157 


0.3156 


0.3163 


In 


163.35 


163.81 


163.87 


163.46 




0.1718 


0.1714 


0.1712 


0.1721 




±3 


±3 


±3 


±4 


k 


0.829 


0.820 


0.819 


0.834 




±7 


±4 


±3 


±11 


Vp 


0.0939 


0.0942 


0.0941 


0.0938 


Ts 


0.0779 


0.0772 


0.0771 


0.0782 


Up 


0.50 


0.59 


0.63 


0.52 


Us 


0.51 


0.59 


0.62 


0.54 


Vp 


0.83 


0.95 


1.07 


1.26 


Vs 


0.84 


0.96 


1.08 


1.28 


Js 


1.0153 


1.0088 


0.9884 


1.0269 




±13 


±12 


±12 


±19 


Ls 1 Lp 


0.6946 


0.6777 


0.6652 


0.7075 


a (mmag) 


5.8 


5.4 


5.5 


8.0 


A^obs used 


818 


815 


820 


820 



Table B.2. Photometric solutions for SWCMa with 
linear limb darkening coefficients adopted from Claret 
(120001) . 



Parameter 


y 


b 


V 


u 


i (°) 


88.34 


88.47 


88.54 


88.32 




±5 


±4 


±4 


±7 


ecosoj 


-0.30295 


-0.30312 


-0.30317 


-0.30311 




±7 


±5 


±5 


±9 


esinu 


±0.0911 


±0.0878 


±0.0870 


±0.0915 




±18 


±15 


±15 


±26 


e 


0.3163 


0.3156 


0.3154 


0.3166 




163.27 


163.84 


164.00 


163.20 


Tp ± Ts 


0.1729 


0.1725 


0.1721 


0.1730 




±3 


±3 


±6 


±5 


k 


0.855 


0.836 


0.830 


0.864 




±12 


±7 


±6 


±19 


Tp 


0.0932 


0.0939 


0.0940 


0.0929 


rs 


0.0797 


0.0786 


0.0781 


0.0802 


Up 


0.60 


0.69 


0.72 


0.59 


Us 


0.61 


0.69 


0.72 


0.60 


yp 


0.83 


0.95 


1.07 


1.26 


ys 


0.84 


0.96 


1.08 


1.28 


Js 


1.0207 


1.0142 


0.9977 


1.0255 




±18 


±18 


±18 


±26 


Ls / Lp 


0.7430 


0.7085 


0.6871 


0.7606 


a (mmag) 


5.8 


5.5 


5.5 


7.9 


A''obs used 


818 


817 


820 


820 



ent LD prescriptions in Table IB. 41 where the results from 
the separate passbands have been averaged in each case, 
with weights inversely proportional to the rms residual of 
each solution. The light elements finally adopted for the 
analysis of SWCMa are those with LD free, and are re- 
peated in Table|3]of Sect. 1^751 in the main text, with more 
conservative errors as described there. 

Appendix C: Light curve solutions for HWCMa 
under different limb-darkening assumptions 

Table IC.ll reports the results from JKTEBOP fits using 
linear LD coefficients from Claret (|2000p . The errors listed 
include the uncertainty in the light ratio constraint, but 
are otherwise internal and unrealistically small in most 
cases. Similarly small errors are obtained using the van 
Hamme p993p coefficients. The geometric elements show 
good agreement between the v, b, and y bands, with u be- 
ing more discrepant (and also more uncertain). Results us- 
ing the LD coefficients from van Hamme (|1993p are given 
in Table IC.2[ and show similarly good agreement in vby. 
In this case we report more conservative errors from 1000 
Monte Carlo simulations in which we perturbed the main 
adjustable quantities that were held fixed to allow their er- 
rors to propagate through: the theoretical LD coefficients 
were allowed to vary by ±0.08, esinw and ecosw were 
perturbed by amounts corresponding to the spectroscopic 
uncertainties in e and w, and the flux ratios Js were al- 
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lowed to vary by ±0.002. Monte Carlo errors using LD 
from Claret (j2000p are very similar to these. 

The systematic differences between the fits with Claret 
and van Hamme LD coefficients are smaller than we found 
before for SW CMa: with the Claret coefficients the incli- 
nation angle is marginally smaller (by 0?01), the sum of 
the radu is ~0.5% larger, and k is also 0.1% larger, all in 
the same direction as found for SW CMa. The individual 
radii are both systematically larger by about 0.5%. 

The light elements we adopt for the analysis of 
HWCMa are those that use the van Hamme (|1993p co- 
efficients. Final values averaged over the vby passbands 
are presented in Table [7| of Sect. 13.31 



Table B.4. Weighted mean photometric solutions for 
SW CMa for three different treatments of the linear limb 
darkening coefficients. 



Parameter 


Free 


van Hamme 


Claret 


* ( ) 


00.09 


88.00 


08.44 




ztlU 


-1-1 n 

ztlU 


e cos Lu 




— U.3U312 


— U.3U3U9 






±y 


-Lin 
±iU 


esincj 


+0.0883 


±0.0889 


±0.0889 




±16 


±15 


±22 


e 


0.3157 


0.3159 


0.3159 




±4 


±4 


±6 


uC) 


163.76 


163.66 


163.65 


±28 


±26 


±39 


Tp + Ts 


0.1718 


0.1715 


0.1726 




±1 


±3 


±4 


k 


0.827 


0.824 


0.843 




±4 


±7 


±15 


Tp 


0.0940 


0.0940 


0.0936 




±1 


±2 


±5 


Ts 


0.0778 


0.0775 


0.0789 




±2 


±5 


±9 



Table C.l. Constrained photometric solutions for 
HW CMa with linear limb darkening coefficients adopted 
from Claret ([20001) . 



Table B.3. Photometric solutions for SWCMa allowing 
the linear limb darkening coefficients to vary freely, with 
the condition that they be the same for the two compo- 
nents. 



Parameter 


y 


b 


V 


u 


i n 


88.53 


88.59 


88.62 


88.62 


±16 


±12 


±11 


±16 


ecosuj 


-0.30300 


-0.30314 


-0.30318 


-0.30321 




±6 


±5 


±5 


±7 


esintj 


±0.0902 


±0.0875 


±0.0869 


±0.0893 




±15 


±14 


±14 


±18 




0.3161 


0.3155 


0.3154 


0.3161 


In 


163.45 


163.89 


164.00 


163.59 


Tp + rs 


0.1720 


0.1719 


0.1717 


0.1716 




±8 


±6 


±6 


±8 


k 


0.832 


0.826 


0.825 


0.824 




±18 


±10 


±8 


±13 


Tp 


0.0939 


0.0941 


0.0941 


0.0941 


rs 


0.0781 


0.0778 


0.0776 


0.0775 


Up — Us 


0.52 


0.64 


0.69 


0.50 




±5 


±4 


±4 


±6 


Vp 


0.83 


0.95 


1.07 


1.26 


Vs 


0.84 


0.96 


1.08 


1.28 


Js 


1.0136 


1.0121 


0.9961 


1.0188 




±36 


±26 


±24 


±50 


Ls 1 Lp 


0.7013 


0.6903 


0.6771 


0.6903 


a (mmag) 


5.8 


5.5 


5.5 


8.0 


A^obs used 


818 


817 


820 


820 



Parameter 


y 


b 


V 


u 


in 


84.79 


84.85 


84.85 


84.66 




±6 


±7 


±5 


±9 


Tp ± rs 


0.06862 


0.06763 


0.06781 


0.07038 




±6 


±8 


±6 


±10 


k 


1.01294 


1.01275 


1.01304 


1.01076 




±4 


±5 


±4 


±5 


Tp 


0.03409 


0.03360 


0.03368 


0.03500 


Ts 


0.03453 


0.03403 


0.03412 


0.03538 


Up 


0.59 


0.67 


0.70 


0.65 


Us 


0.59 


0.67 


0.70 


0.64 


Vp 


0.90 


1.03 


1.17 


1.38 


Vs 


0.88 


1.01 


1.15 


1.36 


Js 


1.075 


1.091 


1.105 


1.078 


Ls/ Lp 


1.103 


1.119 


1.134 


1.106 


a (mmag) 


4.4 


5.5 


4.1 


7.0 


A^obs used 


415 


413 


409 


413 
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Table C.2. Constrained photometric solutions for 
HW CMa with linear limb darkening coefficients adopted 
from van Hamme (|1993|1 . 



Parameter 


y 


b 


V 


u 


i n 


84.80 


84.87 


84.87 


84.67 


±8 


±9 


±8 


±11 


rp + Ts 


0.06829 


0.06725 


0.06741 


0.07008 






±S1 


±fifi 

znuu 


±1 riQ 


If 

K 




1 m 07 

I.UIU ( 


1 ni no 


l.UlUo 




±225 


±223 


±223 


±227 


rp 


0.03393 


0.03345 


0.03352 


0.03485 




±45 


±50 


±49 


±64 




0.03437 


0.03380 


0.03389 


0.03523 




±55 


±57 


±53 


±69 


Up 


0.51 


0.59 


0.62 


0.58 


Us 


0.51 


0.58 


0.61 


0.57 


Vp 


0.90 


1.03 


1.17 


1.38 




0.88 


1.01 


1.15 


1.36 


Js 


1.075 


1.091 


1.105 


1.078 




±2 


±2 


±2 


±2 


Ls / Lp 


1.103 


1.119 


1.134 


1.106 


±20 


±20 


±21 


±20 


a (mmag) 


4.4 


5.5 


4.1 


7.0 


iVobs used 


415 


413 


409 


413 



